含水合物沉积物多相渗流特性研究进展

天然气水合物作为一种储量大、无污染的清洁能源近些年受到了广泛关注. 近20年来,中国进行了较大范围的陆海域天然气水合物储层勘探与储量预测.2017年,中国地质调查局牵头对南海神狐海域的天然气水合物进行了基于降压渗流原理的试验性开采.国内外已进行的水合物试采工程面临着气体产量低、出砂较多等问题,其最主要的原因之一是开发过程中沉积物内复杂多相渗流机理尚不明晰.本文综述了平行毛细管模型、Kozeny模型等广泛应用于天然气水合物开发渗流分析的理论模型,对比分析了水合物开发多尺度渗流过程模拟方法,简述了国内外含水合物沉积物渗透率测试、渗流过程中沉积物物性演变以及水合物开采室内模拟等方面的渗流实验进展,总结了矿场尺度的天然气水合物储层开采过程中产气数值模拟手段,展望了多相渗流模型、储层原位含水合物样品室内测试及结构与物性演化、矿场尺度数值模拟与水平井压裂技术等应用研究的未来方向与挑战.      

Numerical analysis of swirling turbulent droplet-laden flow and heat transfer in a sudden pipe expansion

The effect of swirling intensity on the structure and heat transfer of a turbulent gas–droplet flow after a sudden pipe expansion has been numerically simulated. Air is used as the carrier phase, and water, ethanol, and acetone are used as the dispersed phase. The Eulerian approach is applied to simulate the dynamics and heat transfer in the dispersed phase. The gas phase is described by a system of Reynolds-averaged Navier-Stokes (RANS) equations, taking into account the effect of droplets on mean transport and turbulent characteristics in the carrier phase. Gas phase turbulence is predicted using the second-moment closure. A swirling droplet-laden flow is characterized by an increase in the number of small particles on the pipe axis due to their accumulation in the zone of flow recirculation and the action of the turbulent migration (turbophoresis) force. A rapid dispersion of fine droplets over the pipe cross-section is observed without swirling. With an increase in swirling intensity, a significant reduction in the length of the separation region occurs. The swirling of a two-phase flow with liquid droplets leads to an increase in the level of turbulence for all three types of liquid droplets investigated in this work due to their intensive evaporation. It is shown that the addition of droplets leads to a significant increase in heat transfer in comparison with a single-phase swirling flow. The greatest effect of flow swirling on heat transfer intensification in a two-phase gas-droplet flow is obtained for the droplets of ethanol and water and smallest effect is for the acetone droplets.     

Simulating bubble dynamics in a buoyant system

Multiphase flows are critical components of many physical systems; however, numerical models of multiphase flows with large parameter gradients can be challenging. Here, two different numerical methods, volume of fluid (VOF) and smoothed particle hydrodynamics (SPH), are used to model the buoyant rise of isolated gas bubbles through quiescent fluids for a range of Bond and Reynolds numbers. The VOF is an Eulerian grid–based method, whereas the SPH is Lagrangian and mesh free. Each method has unique strengths and weaknesses, and a comparison of the two approaches as applied to multiphase phenomena has not previously been performed. The VOF and SPH simulations are compared, verified, and validated. Results using two-dimensional VOF and SPH simulations are similar to each other and are able to reproduce numerical benchmarks and experimental results for sufficiently large Morton and Reynolds numbers. It is also shown that at low Reynolds numbers, the two methods, SPH and VOF, diverge in the transient regime of the bubble rise. Regimes that require simulations capable of representing three-dimensional drag are identified as well as regimes in which results from VOF and SPH diverge.     

An improved particle shifting technique for incompressible smoothed particle hydrodynamics methods

The smoothed particle hydrodynamics (SPH) method is one of the powerful Lagrangian tools for modeling free surface flows. However, it suffers from particle disorder, which leads to interpolation and numerical errors. To overcome this problem, several techniques have been introduced until now, among which the particle shifting technique (PST) based on Fick's law is an efficient one. The current form of this method needs tuning parameters to fulfill numerical stability criteria. In this study, to eliminate calibration factors, a new shifting coefficient is derived theoretically based on particle positions before and after shifting, regardless of other parameters such as velocity, pressure, time step intervals, etc. The only required input is particle positions, and the main concern is conserving particle densities in their updated positions. In addition to the proposed PST, a new distribution index (DI) is introduced for measuring the spatial uniformity of particles. Furthering the research, some novel treatments are also studied to improve particle movements near free surface boundary. The proposed idea is only assessed for ISPH method in this study, and its performance in other SPH schemes needs more investigations. Following this innovative method, it is validated by modeling different cases including dam break flow, paddle movement, and elliptical water drop. In all cases, particle arrangements have been improved by means of the modified shifting method. In that sense, good agreements between simulation results with experimental data, analytical solutions, and other numerical methods approve the ability of the developed method in simulating free surface flows.     

A discontinuous Galerkin immersed boundary solver for compressible flows: Adaptive local time stepping for artificial viscosity–based shock-capturing on cut cells

We present a higher-order cut cell immersed boundary method (IBM) for the simulation of high Mach number flows. As a novelty on a cut cell grid, we evaluate an adaptive local time stepping (LTS) scheme in combination with an artificial viscosity–based shock-capturing approach. The cut cell grid is optimized by a nonintrusive cell agglomeration strategy in order to avoid problems with small or ill-shaped cut cells. Our approach is based on a discontinuous Galerkin discretization of the compressible Euler equations, where the immersed boundary is implicitly defined by the zero isocontour of a level set function. In flow configurations with high Mach numbers, a numerical shock-capturing mechanism is crucial in order to prevent unphysical oscillations of the polynomial approximation in the vicinity of shocks. We achieve this by means of a viscous smoothing where the artificial viscosity follows from a modal decay sensor that has been adapted to the IBM. The problem of the severe time step restriction caused by the additional second-order diffusive term and small nonagglomerated cut cells is addressed by using an adaptive LTS algorithm. The robustness, stability, and accuracy of our approach are verified for several common test cases. Moreover, the results show that our approach lowers the computational costs drastically, especially for unsteady IBM problems with complex geometries.     

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

主动射流控制水翼空化的数值模拟与分析

为了改善高速流动工况下水翼吸力面上流场的空化特性,提出了水翼表面主动射流对绕水翼周围流动加以控制的方法.基于密度分域滤波的FBDCM混合湍流模型联合Zwart-Gerber-Belamri空化模型,分析了来流空化数为0.83,来流攻角为8°,射流位置距水翼前缘为x=0.19c时,主动射流对于水翼吸力面上流动的空化特性和水动力特性影响.对回射流的强度进行了量化分析,以探究回射流与流场空化特性的关系.数值分析结果表明,在射流水翼吸力面上的时均空泡体积为原始水翼的1/15,使得流场内空化流动由云空化状态转变为较为稳定的片空化状态,显著地削弱了云空化的发展.此外,射流极大地改善了水翼的水动力性能,使得水翼的升阻比较原始水翼提高了22.9%,空泡的脱落频率减少了26.2%,空泡脱落所引起的振幅减小了9.1%.射流大幅降低了水翼吸力面上低压区面积,水翼吸力面上流体的逆向压力减小,回射流强度降低;同时,射流使水翼吸力面上的边界层减薄,增强了流动的抗逆压梯度能力,一定程度上阻挡了回射流向水翼前缘的流动,这也从机理上分析了主动射流抑制空化的原因.      

船舶与海洋平台结构冰载荷的高性能扩展多面体离散元方法

船舶与海洋平台结构的冰载荷是寒区海洋工程结构物设计中的关键参数,而离散元方法是有效计算结构冰载荷的重要手段. 本文采用基于闵可夫斯基和原理的扩展多面体离散元方法模拟船舶与海洋平台结构的相互作用过程. 其中,构造扩展多面体的近似包络函数并建立了基于优化模型的快速接触搜索算法;考虑单元间粘结作用的刚度软化过程建标识码元间的粘结-破碎模型. 同时,发展了 CPU-GPU 协同异构环境下的高性能并行算法. 为分析海冰与海洋结构作用中的冰载荷,采用ISO标准验证了扩展多面体离散元分析结构冰载荷的准确性. 采用离散元方法计算了船舶结构的冰载荷,研究了船舶结构表明的线载荷分布特点,并采用船舶结构冰阻力经验公式验证了计算结果的合理性. 采用离散元方法计算了平整冰区与多桩腿平台结构的相互作用,分析各桩腿上的冰载荷特点. 针对碎冰区的海冰管理过程,采用离散元方法分析了船舶结构绕行过程中的船舶和海洋平台结构冰载荷. 本文方法可有效应用于海洋结构冰载荷分析,能为极地船舶与海洋平台结构的设计和安全运行提供科学的分析手段.      

THE MODIFIED TRIPLES ALGORITHM FOR FLUSH AIR DATA SENSING SYSTEM APPLIED TO THE VEHICLE WITH BLUNT FORE-BODIES1)

针对传统的用于钝头体飞行器的嵌入式大气数据传感系统的经典三点式算法对测压孔配置约束性严格的特点,对经典三点式算法进行了改进及验证。首先,对于经典三点式算法及改进算法进行了论述;其次,针对不同测压孔配置对算法精度的影响进行了系统的评估;最后,对于传统的经典三点式算法与改进三点式算法的优劣进行了比较。结论为：(1)改进的三点式算法对于测压孔的配置敏感性变差,即对测压孔约束很宽松,可用测压孔数增加,从而使得算法的适用性更强;(2)改进的三点式算法的精度与三点式算法相当,但是需要对测压孔进行系统选取及对比验证。     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界，避免生成复杂贴体网格及动网格，极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法，将简化附加体积力以源项形式嵌入动量方程迭代中，通过用户自定义函数对CFD软件FLUENT二次开发，实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证，并探讨了插值函数对计算精度的影响。研究表明，通过引入浸入边界模型，能够提高计算效率，并实现结构网格背景下复杂边界和动边界的高效建模。